Theory of Ferroelectric Nanoparticles in Nematic Liquid Crystals 



Lena M. Lopatina and Jonathan V. Sehngeiu 
Liquid Crystal Institute, Kent State University, Kent, OH 44^4^ 
(Dated: February 18, 2009) 

Recent experiments have reported that ferroelectric nanoparticles have drastic effects on nematic 
hquid crystals — increasing the isotropic-nematic transition temperature by about 5 K, and greatly 
increasing the sensitivity to applied electric fields. To understand these effects, we develop a theory 
for the statistical mechanics of ferroelectric nanoparticles in liquid crystals. This theory predicts the 
enhancements of liquid-crystal properties, in good agreement with experiments. These predictions 
apply even when electrostatic interactions are partially screened by moderate concentrations of ions. 



In recent years, many experiments have found that col- 
loidal particles in nematic liquid crystals exhibit remark- 
able new types of physical phenomena. If the particles 
are micron-scale, they induce an elastic distortion of the 
liquid-crystal director. This elastic distortion leads to an 
effective interaction between particles, and offers the pos- 
sibility of organizing a periodic array ofparticles, with 
possible photonic applications [l|, |3j B H S Si- If the 
particles are 10-100 nni in diameter, they are too small 
to distort the liquid-crystal director, and hence the sys- 
tem enters another range of behavior. Experiments have 
shown that low concentrations of ferroelectric nanopar- 
ticles can greatly enhance the physical properties of ne- 
matic liquid crystals 0, H, B K [ill H E, [ll • In partic- 
ular, Sn2P2S6 or BaTiOs nanoparticles at low concentra- 
tion (<1%) increase the orientational order parameter of 
the host liquid crystal, and increase the isotropic-nematic 
transition temperature by about 5 K. The nanoparti- 
cles also decrease the switching voltage for the Fred- 
ericksz transition. These experimental results are im- 
portant for fundamental nanoscience, because they show 
that nanoparticles can couple to the orientational order 
of a macroscopic medium. They are also important for 
technological applications, because they provide a new 
opportunity to tune the properties of liquid crystals with- 
out additional chemical synthesis. 

A key question in this field is how to understand and 
control the properties of liquid crystals doped with fer- 
roelectric nanoparticles. One paper has argued that the 
nanoparticles produce large local electric fields, which po- 
larize the liquid-crystal molecules and thereby increase 
the intermolecular interaction [l3j . This increased in- 
teraction leads to a higher isotropic-nematic transition 
temperature Tni- A limitation of this approach is that 
it does not consider the orientational distribution of the 
nanoparticles themselves. One might expect this distri- 
bution to be crucial for understanding how nanoparticles 
affect the orientational order of the liquid crystal. 

In this paper, we propose a new theory for the sta- 
tistical mechanics of ferroelectric nanoparticles in liquid 
crystals, which is based specifically on the orientational 
distribution of the nanoparticle dipole moments. This 
distribution is characterized by an orientational order pa- 
rameter, which interacts with the orientational order of 



the liquid crystals and stabilizes the nematic phase. We 
estimate the coupling strength and calculate the result- 
ing enhancement in Tni, in good agreement with experi- 
ments. This enhancement occurs even when electrostatic 
interactions are partially screened by moderate concen- 
trations of ions in the liquid crystal. In addition, we 
predict the response of the isotropic phase to an applied 
electric field, known as the Kerr effect, and show that it 
is greatly enhanced by the presence of nanoparticles. 

To begin the calculation, consider a spherical nanopar- 
ticle with radius R and electrostatic dipole moment p, 
surrounded by a nematic liquid crystal, as shown in 
Fig. 1. The electric field E generated by the nanopar- 
ticle interacts with the order tensor Qj^^ of the liquid 
crystal through the free energy 
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where Ae is the dielectric anisotropy of the fully aligned 
liquid crystal. The electric field of the nanoparticle has 
the standard dipolar form 
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neglecting higher-order corrections due to the dielectric 
anisotropy of the liquid crystal. Near the nanoparticle, 




FIG. 1: Nanoparticles surrounded by liquid crystal, (a) Par- 
ticle with no electric dipole moment, in the isotropic phase, 
(b) Ferroelectric particle with electric dipole moment, which 
produces an electric field that interacts with the orientational 
order of the liquid crystal. 
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the electric field varies rapidly as a function of position. 
However, the liquid-crystal order cannot follow that rapid 
variation, because it would cost too much elastic energy. 
Hence, for sufficiently small nanoparticles, we can assume 
that order tensor Q]f^ is uniform in space. In that case, 
we can integrate the interaction free energy to obtain 
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Now consider a low concentration pnp of nanoparticles 
dispersed in the liquid crystal. The dipole moments of 
these nanoparticles will not all have the same orientation; 
rather there must be a distribution of orientations. As a 
result, the interaction of liquid crystal and nanoparticles 
gives the free energy density per unit volume 
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averaged over the distribution of nanoparticle orienta- 
tions. This distribution can be expressed in terms of a 
nanoparticle order tensor 
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analogous to the standard liquid-crystal order tensor. 
Hence, the interaction free energy density becomes 
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which shows an explicit coupling between the order ten- 
sor of the liquid crystal and the order tensor of the 
nanoparticles. If we make the reasonable assumption 
that both of these tensors are aligned along the same 
axis, then this interaction reduces to 
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where ^lc and ^np are the scalar order parameters of 
the liquid crystal and the nanoparticles, respectively. 

To model the statistical mechanics of nanoparticles dis- 
persed in the liquid crystal, we must expand the free en- 
ergy in both order parameters S'lc and S'np, which gives 
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Here, the first three terms are the standard Landau-de 
Gennes free energy of a nematic liquid crystal. The first- 
order isotropic-nematic transition of the pure liquid crys- 
tal occurs at Tni = T* + {2b'^) / {9a[^Qc) , and the leading 
coefficient in this expansion can be estimated through 
Maier-Saupe theory as a[^Q = S/cbPlC; where pLC is the 
concentration of liquid-crystal molecules per volume [l5| . 



The fourth term in the free energy is the entropic cost 
of imposing orientational order on the nanoparticles. By 
expanding the entropy in terms of the orientational distri- 
bution function of the nanoparticles, we can estimate the 
coefficient as onp = SfceT/ONP- The final term is the cou- 
pling between the liquid-crystal order and the nanopar- 
ticle order, calculated above. 

We minimize the free energy of Eq. ([8]) over the 
nanoparticle order parameter to find the optimum value 
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This equation shows that the liquid crystal induces ori- 
entational order of the nanoparticles, with a nanoparticle 
order parameter proportional to the liquid-crystal order 
parameter. Note that the induced order is independent 
of the nanoparticle concentration, which is reasonable be- 
cause it arises from interaction of individual nanoparti- 
cles with the liquid crystal, not from interactions between 
nanoparticles. We then substitute this expression back 
into the free energy to obtain 
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In this equation, the coefficient of has been shifted 
by the interaction with the nanoparticles. This shift in- 
creases the isotropic-nematic transition temperature by 
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The last expression has been simplified by writing p — 
{^TrR^)P and pnp = 0Np/(|7ri?'^), where P is the polar- 
ization and (f>^p the volume fraction of the nanoparticles. 

To estimate ATni numerically, we use the following 
parameters appropriate for Sn2P2S6 nanoparticles in the 
liquid crystal 5CB: </>np = 0.5%, i? = 35 nm, Tni = 308 
K, PLC = 2.4 X 10^7 m-3, P = 0.04 Cm'^, fee = 1-38 x 
10-23 JK-i, eo = 8.85 x IQ-^^ C^N-^m-^ and Ae w 
e « 10 [3]. With those parameters, we obtain ATni « 
5 K, which is roughly consistent with the increase that is 
observed experimentally. Of course, there is a substantial 
uncertainty in this estimate, because the parameters R 
and P are not known very precisely in the experiments. 

Note that our model predicts that the enhancement 
ATni should be ffi'st-order in volume fraction ^np, 
fourth-order in polarization P, and third-order in R. In 
particular, increasing R should increase ATni as long 
as the nanoparticles are not large enough to disrupt 
the liquid-crystal order. This prediction disagrees with 
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Ref. [l3|, which predicts that ATni should be first-order 
in second-order in P, and independent of R. 

At this point, we must consider the effects of ionic 
impurities in the liquid crystal. Any liquid crystal con- 
tains some concentration of free positive and negative 
ions, which can redistribute in response to electric fields. 
One might worry that these ions would screen the elec- 
tric field of the nanoparticles, and hence prevent the en- 
hancement of Tni- To address this issue, we solve the 
linearized Poisson-Boltzmann equation to find the elec- 
tric field around a dipole in the presence of ions 
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In this expression, k ^ is the Debye screening length 
given by 
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where n is the concentration and q the charge of the 
ions. With this expression for the field, we can repeat the 
calculation above for the enhancement in Tni, leading to 
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To interpret this result, note that the key parameter 
is kR, the ratio of the nanoparticle radius to the Debye 
screening length. If the ion concentration is low, then the 
screening length is large compared with the nanoparticle 
radius, and hence ATni is as large as in the unscreened 
case. However, if the ion concentration is sufficiently 
high, then the screening length becomes comparable to 
the nanoparticle radius, and hence the enhancement is 
screened away. For a specific example. Fig. 2 shows Tni 
as a function of ion concentration n, using the numerical 
parameters discussed above. In this example, the full un- 
screened enhancement persists up to n « 10^° ions/m^. 
It then decays away as a function of ion concentration, 
and is virtually eliminated by n ^ 10^^ ions/m^. 

Typical measurements of the ion concentration in 5CB 
show n « 10^" ions/m^ and hence 260 nm [l3]. 

Because this screening length is much greater than the 
nanoparticle radius, the enhancement should indeed be 
observable in realistic experiments. Note that the ion 
concentration in a liquid crystal varies over several or- 
ders of magnitude, depending on preparation conditions. 
Hence, we speculate that variations in ion concentration 
may be one explanation for variations in published ex- 
perimental measurements of ATni. 

So far, we have modeled the spontaneous ordering of 
a nanoparticle-doped liquid crystal. We can also use the 
same theoretical approach to predict how the system re- 
sponds to an applied electric field. For a specific example, 
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FIG. 2: Predicted isotropic-nematic transition temperature as 
a function of ion concentration in a nanoparticle-doped liquid 
crystal, using numerical parameters presented after Eq. 



we investigate the Kerr effect, in which an applied electric 
field E induces orientational order in the isotropic phase, 
slightly above the isotropic-nematic transition. In a pure 
liquid crystal, the Kerr effect is a weak alignment propor- 
tional to E^. In a liquid crystal doped with ferroelectric 
nanoparticles, we expect that an applied electric field will 
induce polar order of the nanoparticles, proportional to 
E. This polar order will necessarily induce nematic or- 
der of the nanoparticles, proportional to E^, which will 
in turn induce nematic order of the liquid crystal, also 
proportional to E^. Hence, the nanoparticle-doped liq- 
uid crystal should have an enhanced Kerr effect with the 
same symmetry as the standard Kerr effect, but with a 
much larger magnitude. 

To model the enhanced Kerr effect, we must generalize 
the Landau theory presented above in three ways. First, 
we must introduce a polar order parameter Ma — {pa)/p 
for the nanoparticles, as well as the nematic order param- 
eters Q^J and Q^^. Second, we must consider the ener- 
getic coupling of an applied electric field to the order pa- 
rameters. In the free energy density, an applied field cou- 
ples linearly to the polar order parameter of the nanopar- 
ticles through the interaction —p^ppEaMa, and couples 
quadratically to the nematic order parameter of the liq- 
uid crystal through the interaction —^eoAeEaEpQ^^'j^. 
Third, we must calculate the entropy of a nanoparticle 
distribution characterized by both order parameters Ma 
and Q^p, following the method of Ref. 15]. Assuming 
that all order parameters are aligned along the electric 
field direction, the free energy becomes 



F 
V 



<ciT-T*) , 
^ J' 



LC 



:S\ 



LC 



•-"LC 



eoAe 



' 1807reoe2i?3 
-/sbT/9np 



5'lc5'np - pNPpEM 



E 5lc 



(15) 



5 r.2 

2'-'np 



-M^ - a^NpAT 
2 



4 



0) 
0) 

E 

CO 

i_ 

(0 
Q. 

"D 
6 



1.0 

0.8 

0.6 H 

0.4 

0.2 

0.0 



- V/|im 

- 20 V/|im 

- 40 V/|im 
60 V/nm 
V/nm 
0.02 V/^im 
0.04 V/^im 
0.06 V/^im 



290 300 310 320 330 
Temperature (K) 



340 350 



FIG. 3: (Color online) Prediction for field-induced order pa- 
rameter Slc as a function of temperature for several values of 
applied electric field, with and without ferroelectric nanopar- 
ticles, using numerical parameters presented after Eq. 



nanoparticle doping at and above the isotropic-nematic 
transition, in the regime where Landau theory is vaUd. 

In conclusion, we have developed a theory for the sta- 
tistical mechanics of ferroelectric nanoparticles in ne- 
matic liquid crystals. This theory predicts the enhance- 
ment in the isotropic-nematic transition temperature and 
in the response to an applied electric field, which can 
be tested experimentally. The work demonstrates the 
coupling of nanoparticles with macroscopic orientational 
order, and provides an opportunity to improve the prop- 
erties of liquid crystals without chemical synthesis. 

We would like to thank J. L. West, Y. Reznikov, and 
P. Bos for many helpful discussions. This work was sup- 
ported by NSF Grant DMR-0605889. 



We minimize this free energy over all three order pa- 
rameters, M, 5np, and S'lc- In the high-temperature 
isotropic phase, in the limit of small electric field, the 
resulting liquid-crystal order parameter is 
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In this expression, note that the induced order parameter 
depends on electric field and temperature exactly as in 
the standard Kerr effect, but the coefficient is increased 
by the coupling with nanoparticles. In the square brack- 
ets, the first term of 1 indicates the standard Kerr effect 
for pure liquid crystals, and the second term indicates 
the relative enhancement due to nanoparticle doping. 

For a specific numerical example, we use the same pa- 
rameters presented after Eq. (fTTjl . With these parame- 
ters, the relative enhancement in the Kerr effect is ex- 
tremely large, of order 10''. Figure 3 plots the predicted 
order parameter S'lc a function of temperature for sev- 
eral values of the applied electric field, with and without 
nanoparticles. This plot shows explicitly that the pres- 
ence of nanoparticles greatly enhances the sensitivity to 
applied electric fields in the isotropic phase, as well as 
enhancing the isotropic-nematic transition temperature. 
This prediction should be tested in future experiments, 
and should provide an opportunity to build liquid-crystal 
devices that can operate at lower electric fields. 

As a final point, we should mention one limitation of 
our model. Like all Landau theories, our model involves 
an expansion of the free energy in powers of the order 
parameters, and hence it overestimates the order param- 
eters that occur in the low-temperature phase. Future 
work may extend this model through asymptotic low- 
temperature approximations to the free energy. Never- 
theless, our current model clearly shows that effects of 
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